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Abstract 

We measure the complex scattering amplitudes of a flat microwave cavity (a "chaotic billiard"). 
Time-reversal (7~) invariance is partially broken by a magnetized ferrite placed within the cavity. 
We extend the random-matrix approach to T violation in scattering, determine the parameters 
from some properties of the scattering amplitudes, and then successfully predict others. Our work 
constitutes the most precise test of the theoretical approach to T violation within the framework 
of random-matrix theory so far available. 

PACS numbers: 24.60.Ky, 05.45.Mt, ll.30.Er, 85.70.Ge 
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We measure the effect of partial violation of time-reversal (T) invariance on the exci- 
tation functions of a flat microwave cavity induced by a magnetized ferrite placed within 
the cavity. The classical dynamics of a point particle moving within the cavity and elas- 
tically reflected by the walls, is chaotic. The statistical properties of the eigenvalues and 
eigenfunctions of the analogous quantum system are, therefore, expected to follow random- 
matrix predictions Random-matrix theory (RMT) provides a universal description of 
generic properties of chaotic quantum systems. In particular, RMT yields analytical expres- 
sions for correlation functions of scattering amplitudes [2j that can be generalized to include 
Tviolation. Although widely used (to discover signatures of Tviolation in compound-nucleus 
reactions Q] in the Ericson regime \&\, to describe electron transport through mesoscopic 

n 

samples in the presence of a magnetic field [51], and in ultrasound transmission in rotational 
flows [6]), that generic model for Tviolation has, to the best of our knowledge, never been 
exposed to a detailed experimental test. With our data we perform such a test. 

Our aim is not a detailed dynamical modeling of the properties of the cavity. With the 
exception of the average level density we determine the parameters of the RMT expressions 
from fits to some of the data. We then test the RMT approach by using it to predict 
other data, and by subjecting our fits to a thorough statistical test. All of this is in the 
spirit of a generic RMT approach since a dynamical calculation of the relevant parameters 
is not possible for many systems. Such a calculation works only for special chaotic quantum 
systems like some cavities where the semiclassical approximation can be used [3, Q] ■ We use 
that approximation only to determine the average level density, and to estimate the range 
of validity of RMT in terms of the shortest periodic orbit. 

?een used before to study the effect of T-invariance violation 
11 ] and on the eigenfunctions jl2| . Here we study fluctuations of 
the scattering amplitudes versus microwave frequency. For our cavity the average resonance 
spacing d is of the order of the resonance width T, and we work in the regime of weakly 
overlapping resonances. 

n 

Experiment. The flat copper microwave resonator has the shape of a tilted stadium 131 ] 
(see Fig. [I]) and a height of 5 mm. The excitation frequency / ranges from 1 to 25 GHz. In 
that range, only one vertical mode of the electric field strength is excited. The Helmholtz 
equation for the tilted stadium is then mathematically equivalent to the Schrodinger equation 
of a two-dimensional chaotic quantum billiard 14j. An Agilent PNA-L N5230A vector 
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network analyzer (VNA) coupled rf energy via one of two antennas labeled 1 and 2 into the 
resonator and determined magnitude and phase of the transmitted (reflected) signal at the 
other (same) antenna in relation to the input signal and, thus, the elements S a b(f) with a,b = 
1, 2 of the complex- valued 2x2 scattering matrix S(f). Distorting effects of the connecting 
coaxial cables were removed by calibration. We measured the elements of S(f) in the 
frequency range 1-25 GHz at a resolution of 100 kHz. To improve the statistical significance 
of the data set, an additional scatterer (an iron disc of 20 mm diameter) was placed within 
the cavity. It could be freely moved and allowed the measurement of statistically independent 
spectra, so-called "realizations". 

Time-reversal invariance is violated [ijj] by a ferrite cylinder (Air Ms = 1859 Oe, AH = 
17.5 Oe, courtesy of AFT Materials GmbH, Backnang, Germany) of 4 mm diameter and 
5 mm height. The cylinder was placed inside the resonator and magnetized by an external 
magnetic field B. The field was provided by two NdFeB magnets (cylindrical shape, 20 mm 
diameter and 10 mm height) attached from the outside to the billiard. Field strengths of 
up to 360 mT could be attained. Here we focus on the results at B — 190 mT as there the 
effects are most clearly visible. The spins within the ferrite precess collectively with their 
Larmor frequency about the external field. The rf magnetic fields of the resonator modes 
are, in general, elliptically polarized and couple to the spins of the ferrite. The coupling 
depends on the rotational direction of the rf field. An interchange of input and output 
channels changes the rotational direction and thus the coupling of the resonator modes to 
the ferrite. Figure [2] demonstrates that reciprocity, defined by Su(f) = 52i(/) and implied 
by T invariance, is violated. 

As a measure of the strength of T-invariance violation, we define the cross-correlation 
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FIG. 1: The tilted stadium billiard (schematic). The two antennas 1, 2 connect the resonator to 
the VNA. The ferrite is fixed, the scatterer can be moved freely. 
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coefficient C, 



() ,, ss (e = 0) where 

Cc3 



(1) 



V(\Sn(f)\ 2 )(\S 21 (fW)' 

If T invariance holds, we have C C ross(0) — 1 while for complete breaking of T invariance S\2 
and S21 are uncorrelated and thus C croS s(0) = 0. The average ( ■ ) over the data is taken in 
frequency windows of width 1 GHz and over 6 realizations, i.e., positions of the additional 
scatterer. The upper panel of Figure [3] shows C C ross(0) for the different frequency windows. 
The cross-correlation coefficient is seen to depend strongly on / although complete violation 
of T invariance is never attained. At 5-7 GHz the Larmor frequency of the ferrite matches 
the rf frequency, and the ferromagnetic resonance directly results in C crO ss(0) ~ 0.8. Around 
15 GHz the effects of T-invariance violation are strongest, C cross (0) ~ 0.5. A third minimum 
is observed at about 24 GHz. The connection of the latter two minima to the properties of 
the ferrite is not clear. 

Analysis. We analyze the data with a scattering approach developed in the context of 
compound- nucleus reactions [if]]. The scattering matrix for the scattering from antenna b 



to antenna a with a, b = 1, 2 is written as 



SM) = ^ab ~ 2vr i (H/t (/ - H^y 1 W) t 



(2) 



The matrix W^ a is rectangular and describes the coupling of the N resonant states \i in the 
cavity with the antennas a = 1,2. We assume that T-invariance violation is due to the ferrite 
only. Then W^ a is real. The resonances in the cavity are modeled by H cS = H — iirWW^. 
Here H is the Hamiltonian of the closed resonator. The elements of the real matrix W^ c are 
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FIG. 2: Transmission spectra for B = 190 mT in the range 16-17 GHz. The amplitudes and phases 
of S\2 (solid) and S21 (dashed) are seen to differ. 



4 




5 10 15 20 

Frequency (GHz) 



25 



FIG. 3: Experimentally determined values of C croS s(0) (upper panel) from Eq. (JTJ) and the parameter 
£ for T-invariance violation deduced from these (lower panel) with the help of Eq. ([3]). The error 
bars indicate the r.m.s. variation of C croS s(0) over the 6 realizations. 

equal to those of for c = 1,2. As done successfully before [l?], Q, Ohmic absorption 



of the microwaves in the walls of the cavity and the ferrite is mimicked [19] by additional 
fictitious weakly coupled channels c. The classical dynamics of a point particle within the 
tilted stadium billiard is chaotic. Therefore [l], we model H by an ensemble of random 



matrices. The iV-dimensiona 
as the sum of two parts 



Hamiltonian matrix H of the system (the cavity) is written 
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221], H = H s + i(-K£/VN)H a . The real, symmetric, and 



T-invariant matrix H s is taken from the Gaussian Orthogonal Ensemble (GOE) while the 
real, antisymmetric matrix H a with Gaussian-distributed matrix elements models the T- 
invariance breaking part of H. For 7r£ / yN = 1 the Hamiltonian H belongs to the Gaussian 
Unitary Ensemble (GUE) describing systems with complete T breaking. However, for iV — > 
oo, Tinvariance is significantly broken already when the dimensionless parameter £ is close 
to unity In the same limit C croS s(0) in Eq. (OQ) can be expressed analytically in terms of 
a threefold integral involving the parameter £. For the derivation we extended the method 
of Ref. (24]] where the ensemble average of \S a b\ 2 was computed as function of the parameter 
£. The cross-correlation coefficient C cross (0) is obtained by setting e = 0, o = —1, and 
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a, b = 1, 2 in the function 
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with the notations 



(3) 
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ft 



7T £ , n = A(x + X 1 )(x + X 2 ), 

2^X1(1 + xi)x 2 (l + x 2 ), T = 4x(l - x) 

l±exp(-2tF), A = l-2x, 

V(l + xi) (1 + x 2 ) + xix 2 - (-1)* W, 



A 



l,i = 1,2. 



(4) 



The integration measure fj,(x, Xi, x 2 ) and the function J a b(x,xi,x 2 ) are given explicitly in 



Ref. 



121 



while the 



'unctions K ab (X, X u A 2 |T a , T b , £) and L ab (X, X u X 2 \T a , T b , £) can be read off 
Eq. (2) of Ref. |25j. We checked our analytic results by numerical RMT simulations. The 
parameters of Eq. for e = are £, the transmission coefficients T a = 1 — ((S^a)] 2 for a 



1,2, and the sum r a b s of 300 transmission coefficients that model the Ohmic losses [18|, 1191 ] . 

For a typical set T 1; T 2 , r abs , Fig.@]shows C cross (0) versus £. Within the frequency range 1- 
25 GHz, C croS s(0) depends very weakly on Ti, T 2 , r a b s , and Fig.|4]can be taken to be universal. 
For each data point shown in the upper panel of Fig. [3] the corresponding value of £ was 
read off Fig. H] and the result is shown as function of / in the lower panel of Fig. [3J Figure [3] 
shows that in the interval from 1 to 25 GHz, the ratio of the average resonance width Y to 
the average resonance spacing d varies from T/d ~ 0.01 to Y/d ~ 1.2 while the strength 
£ of T breaking varies from zero to 0.3. Numerical calculations show that for £ = 0.3 the 
spectral fluctuations of the Hamiltonian H for the closed resonator defined below Eq. ([2j) 
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almost coincide with those of the GUE [26fl. We also found that for £ = 0.4 they do not 

n 

differ significantly from those presented in Ref. [9j, where the conclusion was drawn, that 
complete T breaking is achieved. However, even for £ = 0.4 the value of C cross (0) is still far 
from zero. This shows that C cross (0) is a particularly suitable measure of the strength £ of 
T violation. 

Autocorrelation function. Since C croS s(0) depends only weakly on the values of T%, T2 and 
Tabs, we used the autocorrelation function C ab (e) for a more precise determination of these 
parameters, especially of r abs . The function 

C ab (e) = (S ab (f) S* ab (f + c)> - |(^(/))| 2 (5) 



was calculated analytically with the method of Ref. |24j as a function of Ti, T 2 , r a b s , £, and d 
and is obtained from Eq. by setting a = +1. It interpolates between the wel 



•known re- 
sults for orthogonal symmetry |2| (full T invariance) and for unitary symmetry 27| (complete 



violation of T invariance). The mean level spacing d was computed from the Weyl formula. 
The Fourier transform of the function C ab (e) was then fitted to the data as in Ref. Q. As 
starting points we used the values of Ti and T 2 obtained from the measured values of S aa (f) 
and of £ determined from C crO ss(0). For each of the 6 realizations the spectra of S ab (f) were 
divided into intervals Af of 1 GHz length. In each interval the Fourier transform C ab (tk) 
of the autocorrelation function ([S]) was calculated for values of tk between 5 ns and 200 ns. 
The lower limit is determined by the length of the shortest periodic orbit in the classical 
billiard; for smaller values of tk the Fourier coefficients are nongeneric [§]. At tk ~ 200 ns the 
values of C ab {t k ) have decayed over more than three orders of magnitude, and noise limits 
the analysis. The time resolution was 1/A/ = 1 ns. We measured four excitation functions 
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FIG. 4: Dependence of the cross-correlation coefficient C croS s(0) on the parameter £ as predicted 
by the random-matrix model for partial violation of T invariance. Also shown is how C cross (0) = 
0.49(3) translates into £ = 0.29(2). 
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Sab{f) taking (a, b) = (1, 1), (1, 2), (2, 1), (2, 2 yielding a total of 4800 Fourier coefficients 
for each interval. For / > 10 GHz the fitted values for T\ and T2 differ by not more than 
7 % from the initial ones. (For smaller / the intervals of 1 GHz width comprise only few 
resonances). The spread of the data is large, see the left panel of Fig. Going to the time 
domain is useful since the S a b(f) are correlated for neighboring / whereas the correlations 
are removed in the ratios of the experimental and the fitted values for Cabifk)- The latter 
are stationary and fluctuate about unity. Thus the statistical analysis is much simplified. 
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FIG. 5: Autocorrelation function for S\2 i n the range of 16-17 GHz and at B = 190 mT. In the 
time domain (a) the data (dots) scatter around the theoretical fit (solid) for T\ = 0.37, T2 = 0.41, 
r abs = 2.9 and £ = 0.25. Transforming the results back into frequency domain (b) confirms the 
good agreement between data and theory. We observe that neighboring data points in (b) are 
correlated, whereas those in (a) are not. 

For each realization the parameters r a b s and £ were obtained by fitting the analytical 
expression for C a b(tk) to the experimental results. The values of £ determined from these 
fits agree with the ones found from the cross-correlation coefficient. To reduce the spread 
we combined the data from all realizations within a fixed frequency interval. The result 
was analyzed with the help of a goodness-of-fit (GOF) test (see Ref. |l8|) that distinguishes 
between full, partial, and no violation of Tinvariance. We defined a confidence limit such 
that the GOF test erroneously rejects a valid theoretical description of the data with a 
probability of 10 %. With this confidence limit the test rejects the fitted expressions for 
Cab{tk) m only 1 out of the 24 available frequency windows or in 4.2 % of the tests. Thus, 
the RMT model correctly describes the fluctuations of the S'-matrix for partial violation of 
Tinvariance in the regimes of isolated and weakly overlapping resonances. 

Elastic enhancement factor. As a second test of the theory we use the values of £ ob- 
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tained from the cross-correlation coefficients (see Fig. [3]) and the parameters T a , TJ,, r a b s 
resulting from the fit of C a b{tk) to predict the values of the elastic enhancement factor 



W 



1/2 



(l^nl 2 ) (l^l 2 ) / (l^l 2 ) with S% = S ah - (S ah ) as a function of /. We use that 



W = a/Cii(0) C22(0)/Ci2(0), see Eq. (J5]). For T-invariant systems, the elastic enhance- 
ment factor decreases from W = 3 for isolated resonances with many weakly coupled open 
channels to W = 2 for strongly overlapping resonances (r > ti). The corresponding values 
for complete violation of Tinvariance are W = 2 and W = 1, respectively 28j. Figure [6] 
compares the analytic results for the enhancement factor W (filled circles) to the data (open 
circles). For small / (where r/d C 1 and £ ~ 0) the experimental results differ from the 
prediction W = 3. Here only few resonances contribute and the errors of the experimental 
values for W are large. Moreover W is determined from only a single value C a b(0) of the 
measured autocorrelation function while the analytic result is based on a fit of the complete 
autocorrelation function. As / increases so does T/d, and W takes values well below 3. 
At frequencies where £ is largest W drops below the value 2 as predicted, a situation that 
cannot arise for T-invariant systems. The overall agreement between both data sets above 
rs 10 GHz corroborates the confidence in the values of £ deduced from the cross-correlation 
coefficients. 
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FIG. 6: Comparison of elastic enhancement factors. For the evaluation of W the autocorrela- 
tion coefficients C a b(0) were determined either directly from the data (open circles) or from the 
analytic result for partial violation of T invariance (filled circles) with £ determined from the cross- 
correlation coefficient. The error bars indicate the variations within the 6 realizations. The dashed 
horizontal lines indicate the limits of W for T-invariant systems in the regime of isolated (W = 3) 
or overlapping (W = 2) resonances. 



Summary. We have investigated partial violation of Tinvariance with the help of a 
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magnetized ferrite placed inside a flat microwave resonator (a chaotic billiard) with two 
antennas. We measured reflection and transmission amplitudes in the regime of isolated 
and weakly overlapping resonances in the frequency range from 1 to 25 GHz and determined 
the cross-correlation function, the autocorrelation functions, and the elastic enhancement 
factor from the data. The results were used as a test of random-matrix theory for scattering 
processes with partial T violation. That theory yields analytic expressions for all three 
observables. The parameters of the theory (Ti,T 2 ,r a b s and the parameter £ for T violation) 
were partly obtained directly from the data but improved values resulted from fits to the 
autocorrelation function. We find that < £ < 0.3. The validity of the theory was tested 
in two ways, (i) A goodness-of-fit test of the Fourier coefficients of the scattering matrix in 
frequency intervals of 1 GHz width yielded excellent agreement, (ii) The elastic enhancement 
factor predicted from the fitted values of the parameters shows overall agreement with the 
data for frequencies above 10 GHz where the experimental errors are small. We conclude that 
the random-matrix description of S'-matrix fluctuations with partially broken Tinvariance 
is in excellent agreement with the data. 
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